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Abstract 

We analyze the self-averaging properties of time-reversed solutions of the paraxial wave equa- 
tion with random coefficients, which we take to be Markovian in the direction of propagation. 
This allows us to construct an approximate martingale for the phase space Wigner transform of 
two wave fields. Using a priori L^-bounds available in the time-reversal setting, we prove that the 
Wigner transform in the high frequency limit converges in probability to its deterministic limit, 
which is the solution of a transport equation. 

1 Introduction 

In time-reversal experiments a signal emitted by a localized source is recorded at an array of trans- 
ducers. It is then re-emitted into the medium reversed in time, that is, the part of the signal that 
is recorded first is sent back last. Because of the time-reversibility of the wave equation the back- 
propagated signal refocuses approximately at the location of the original source because the array 
is limited in size. A striking experimental observation is that the presence of inhomogeneities in 
the medium improves the refocusing resolution. The explanation for this super-resolution is multi- 
pathing: waves in complex media that are captured by the recording array have undergone multiple 
scattering making it effectively larger than its physical size. Another important feature of super- 
resolution in time reversal is that the refocused signal does not depend on the realization of the 
random medium. That is, the refocused signal is deterministic. Super-resolution and self-averaging 
of refocused signals in complicated media has been observed both in laboratory experiments (see 
reviews |l^ and references therein) and in underwater acoustic wave propagation over long dis- 
tances (tens of kilometers) p^ , |20[ |. Time-reversal techniques have numerous applications ranging 
from medicine to communications and, more recently, imaging in random media [Q, 10, |l^. 



The first mathematical analysis of time-reversal in random media was given by Clouet and Fouque 
[ 11 1 , who analyzed refocusing and self-averaging of time- reversed pulses in a one-dimensional layered 



random medium. Their result was extended to a three-dimensional layered medium in |14|. Super- 



resolution in spatial refocusing and its statistical stability for multi-dimensional waves in random 



media was analyzed in [^, 22], in a remote-sensing regime where the paraxial or parabolic wave 
equation can be used. The refocusing of the average signal in a full three-dimensional medium, in 
the regimes of random geometrical optics and radiative transfer (transport), was studied in H, ^. 
We also mention that another source of multipathing is the mixing of waves by the boundaries in 



an ergodic cavity. This has been studied experimentally in |1£] and mathematically in |6| 
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The purpose of this paper is to analyze time reversal in the radiative transfer regime using the 
parabolic wave equation, when the waves interact fully with the random inhomogeneities. We prove 
mathematically that the refocused signal is self-averaging, which means that it does not depend 
on the realization of the random medium. The mathematical quantitiy that we analyze is the 
Wigner measure of a pair of oscillatory solutions of the random Schrodinger equation. In the present 
setting, the random potential depends in a Markovian way on the variable 2, the main direction 
of propagation of the waves. This allowed us to use in Q a martingale method to prove that the 
average of the Wigner distribution converges to a solution of the radiative transfer equation. In 
this paper we use additional regularity of the Wigner measure, available in time-reversal when there 
is some blurring at the recording array, to show that the whole Wigner distribution, and not only 
its average, converges weakly, as a Schwartz distribution and in probability, to the deterministic 
solution of the transport equation. The blurring at the recording array provides a priori bounds for 
the Wigner transform in L^. These bounds and the Markovianity of the random potential in the 
direction of propagation make the time-reversal problem more tractable mathematically and allow 
us to prove in a fairly simple and straightforward manner self-averaging of the time-reversed signal. 

We recall that the Wigner transform is a convenient tool to analyze high frequency wave prop- 
agation in deterministic [^, 21, 25] and random media [25|. Introduced by Wigner in 1932 |29], 
it has been used extensively in the mathematical literature recently. Convergence of the average 
Wigner distribution to the solution of the radiative transfer equation was first proved by H. Spohn 
in |26| for time-independent potentials on small time intervals. This result was extended to global in 
time convergence by L. Erdos and H.-T. Yau [|l^ ]. These proofs involve infinite Neumann (diagram- 
matic) expansions for the solution of the Schrodinger equation and are quite involved technically. 
The corresponding problem with time-dependent potentials is much simpler mathematically. It was 
treated by us in Q in the Markovian case, and by F. Poupaud and A. Vasseur |23| in the case 
of finite-range time correlations. In this paper we use the fact that the Wigner family arising in 
time-reversal problems is more regular than the usual one because blurring is added at the recording 
array. This provides some additional regularity usually obtained by considering mixtures of states 
as, for instance, in |21, 23, 26]. 

The paper is organized as follows: we describe the scaling and obtain an expression for the back- 
propagated signal in terms of the Wigner transform in Section The main result and assumptions 
on the random medium are formulated in Section ^. The proofs are presented in Section |^. 

This work was partially supported by ONR grants N00014-02- 1-0088 and N00014-02-1-0089. 
G. Bal was supported in part by NSF grant DMS-0072008, G. Papanicolaou by grants AFOSR 
F49620-01- 1-0465 and NSF DMS-9971972, and L. Ryzhik in part by NSF grant DMS-9971742 and 
the Alfred P. Sloan Foundation. 



2 The back-propagated signal in the parabohc approximation 

2.1 The back-propagated signal 

The pressure field p{z, x, t) satisfies the scalar wave equation 

1 d'^p 



c^{z, x) dt"^ 



Ap = 0. (1) 



Here c{z, x) is the local wave speed that we will assume to be random, and the Laplacian operator 
includes both direction of propagation, z, and the transverse variable x S M*^. In the physical 
setting, we have d = 2. We consider dimensions d > 1 to stress that our analysis is independent 
of the dimension. If we assume that at time t = 0, the wave field has a "beam- like" structure in 
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the z direction, and if back-scattering may be neglected, we can replace the wave equation by its 



parabolic approximation |27]. More precisely, the pressure p may be approximated as 



where ijj satisfies the Schrodinger equation 

9V 



(2) 



2iK;^^(z,x, k) + AxV'(2,x, k) + K^(n^(z,x) — l)'0(z,x, k) = 0, 
ip{z = 0,x, k) = ipo(x, k) 



(3) 



with Ax the transverse Laplacian in the variable x. The refraction index ?i(z,x) = co/c(z,x), and 
Co in @ is a reference speed. The rigorous passage to the parabolic approximation from the wave 
equation has been analyzed in |5[ in the deterministic case, and |Q] in a one dimensional random 
medium. A formal derivation of the paraxial approximation that leads to the radiative transfer 
regime is given below in section p.2| . 

We assume that the original source is located in the plane z = 0, and the time-reversal mirror is 
located in the plane z = L as depicted in Figure 




Figure 1: Geometry of the time-reversal experiment. 



During the first stage of a time reversal experiment the signal propagates for a time T, or 
equivalently, over a distance L = cqT, so that the signal arriving at the time-reversal mirror is given 
by 



'0_(L,x)=/ G(L,x,K;y)'0o(y)c?y, 



where the Green's function solves 
dG 

'^iK—{z,x,K;y) AxG(z, x, k; y) +K^{n^{z,x.) - l)G(z, x, y) = 0, 
G(0,i,K;y) =5(x-y). 

Then the signal is time reversed. For three-dimensional acoustic pulses, this means that the pressure 
field is kept unchanged and that the sign of its time derivative is reversed. In the parabolic approx- 
imation, this is equivalent to phase conjugation, where is replaced by its complex conjugate ip*. 
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We assume that the recording array occupies a compact subset of the plane z = L, and introduce 
a real- valued aperture function x{^)- It represents the restriction of the signal onto the array, and 
possible amplification by the array that may vary from one receiver to another. In the absence of 
amplification it is given by the characteristic function of an array set J7 C M*^. We also allow for 
some blurring of the recorded signal, modeled by a convolution with kernel /(x). The signal at L 
after time reversal takes then the form 

jp+{L,x,K)= / x(x)/(x - y)7/;l(L,y;K)x(y)(iy 

= [ x(x)G*(L,y,K;y')x(y)/(x-y)^S(y',^)rfy^y'- 

The last step consists in letting the signal propagate back to the origin z = and time reversing it 
one more time 

V'^(x,k)= / G*{L,x,K;'n)G{L,y,K;y')x{v)x{y)f{v-y)My','^)dydy'dri. (5) 

The last step, phase-conjugating at the source, is not performed in real physical experiments but is 
convenient if we need to compare the back-propagated signal to the original signal. It affects neither 
the degree of refocusing nor the self-averaging effect. The back-propagated signal in time is given 
at the plane z = by 

V;-f^(x,t)= [ e-^'=«''V^(x,K)dK (6) 

JR 

= [ G*{L, X, k; T7)e-^^o^*G(L, y, k; y')x{v)x{y)f{v " y)V'o(y', K)dydy' diqdK. 

We will be interested in the sequel in the refocusing and self- averaging properties of the back- 
propagated signal ip^{'K,K,) for each fixed frequency. Self-averaging of the time signal Tp^{'x.,t) then 
follows from (|6|) when ijjoiy, k) has compact support in k. An interesting conclusion of this paper is 
that the additional blurring by convolution with the kernel / makes the signal ip^{'K, k) self-averaging 
for every frequency k. This should be contrasted with the situation studied in Q, where no blurring 
was introduced and self-averaging was observed only for the full time signal. 



2.2 Localized Source and Scaling 

We recall the formal passage from the reduced wave equation to the parabolic equation (^) described 
in the Appendix to |Q] and explain how the radiative transfer scaling arises in this context. We start 
with the reduced wave equation 

Ap + K^n^{z,x)p = 0, (7) 
and look for solutions of (0) in the form p{z,x.) = e^'^^ip{z,:>c). Then we obtain 

^ + 2iK^ + A^iP + K^n'^ - l)iP = 0. (8) 
oz^ oz 

The refraction index x) is weakly fluctuating so that it has the form 



n2(z,x) = 1 - 2aV 



z x 

r' r 
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where V models random fluctuations that will be described in detail in section Here, Ix and Iz 
are the correlation lengths of V in the transverse and longitudinal directions, respectively, and the 
small parameter a measures the strength of the fluctuation. The waves propagate over a distance 
Lx in the x-plane and in the z-direction, and we rescale x and z accordingly. We also introduce a 
carrier wave number kq and replace k = kqk', k! being the non-dimensional wave number (we drop 
the prime below). The physical parameters determined by the medium are the length scales Ix, Iz 
and the non-dimensional parameter cr ^ 1. 

We now explain the scaling of the parameters Lx, Lz and kq in the radiative transfer regime. 
Equation (jsl) in the non-dimensional variables z' = z/Lz, x' = x/L^^ becomes after we drop the 



primes 



1 d'^ljj 2iKKo dip I 2 2 ( ^Lx X I . 

H — ^ + y^Ax'0-2K^K5cry '0 = 0. (9) 



Li dz'^ Lz dz L| 
We introduce the following small parameters: 



Ix = T' 



Lix 1^0'' X I^Q'' 



and rewrite (0) as 



X 



Let us assume now that 

„2 



5x = 5z<€.l, 72 = 7^ < 1, a = -iz\/ 4 (11) 
and denote e = bx- Then ( |lO| ) becomes after multiplication by e/2 



7..^ d-^ ^ ^J_^ ^ e^^^^ _ p ^ XX ^ ^^^^ 



2 dz 2 ^ Ve' e 

Observe that when k = 0(1) and 7^ <C 1, the first term in (^) is small and may be neglected in 
the leading order since {e'^ipzzl = 0(1). Then (p^) becomes 



iKE^ + ^A,V - ^^V^V (^,-)i; = (13) 



dz 2 \e eJ 

which is the parabolic wave equation (^) in the radiative transfer scaling. 



The second condition in (|r^) implies that the carrier wave number should be chosen as kq = lz fix- 
Then 7^ = Ix/l'i ^ 1 implies that Ix <C lz- Therefore the correlation length in the longitudinal 
direction z should be much longer than in the transverse plane x. Furthermore, we should have 
e = {a l^z)^ = a'^l'^/lx <C 1 which in turn implies that the fluctuation strength a <C (Ix/lz)'^- Given 



that these constraints are satisfied the last condition in (11) implies that the spatial scales Lx and 
Lz should be chosen according to 

T — 1 ^ T — 1 

^X — ''X 9|/L ) ^Z — '•Z 



a 



Note that the constraint lx <^ h implies that Lx <^ Lz, which is the usual constraint for the validity 
of the parabolic approximation. The connection between the relations of the physical parameters 



and the resulting scaling of the parabolic approximation is further discussed in [22|. The rigorous 
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passage to the parabolic approximation in one dimension in a random medium in a similar scaling 
is discussed in 

To define what we mean by the quality of the refocused signal, we assume that the initial pulse 
is centered around a point xq with support comparable to the transverse correlation length. Then 
the initial data for the Schrodinger equation (13) becomes in the rescaled variables 

V;(z = 0,x,k) = Vo(^^^,'^)- (14) 

The transducers should obviously be capable of capturing signals of frequency so that the blurring 
occurs on the scale of the source. We therefore replace /(x) by e~'^/(x/e). Finally, we are interested 
in the refocusing properties of 'tp^{'x.) in the vicinity of xq and consequently introduce the scaling 
X = Xq + e^. We then recast (^) as 

^f(^,/«;xo)= / G*(L,xo + e^,K;?7)G(L,y, k;xo +ey')x('?,y)^o(y',K)dy'(iyd77, (15) 



where 



x{v,y) = x{v)x{y)f{^^). (i6) 



We observe that G(L,x, K;y) = G(L,y,K;x), so that 

V'f(^,K;xo)= / G*(L,xo + e^,K;T7)G(L,xo + ey',K;y)x(^?,y)^o(y','«)dy'dydr7. (17) 

Let us now introduce the auxiliary function as in ^, 

Qe(i,x,Ac;q) = / G{L,^,K;y)x{y)e-'^-''/'dy, (18) 

which satisfies the initial value problem 

i£K-^{z,x; k) + yAxQe(2;,x; k) - K^^/eV (^-, -j Qs{z,:>c,k) = 0, (19) 

Q,(z = 0,x;K)=x(x)e"*''-/^ 

It physically describes the component with wave vector q that is sent back to the plane ^ = by 
the array. We then define the Wigner transform as 

W,iL,x,k,K)= [ /(q)?7,(L,x,k,K;q)dq, (20) 

where 

C/,(L,x,k,K;q) = / e^l'-yQ,(L, x - ^, ac; q)Q*(L, x + ^, a^; q)-^. (21) 

The Wigner transform We is more regular than Ue because of additional averaging in the wave vector 
q. More, precisely, while the family Ue is not uniformly bounded in L^, the family We has a uniform 
bound in L^. This regularizing effect is essentially the same as the one obtained by considering 
mixtures of states for the Schrodinger equation as in |21, 23, p^ ]. 
We can then recast ([T7|) as 



(^,k;xo)= / e^^<^-y)We{L,xo + e^,k,K)My,^)^- (22) 



The above formula shows that the asymptotic behavior of V'e (^j xq) as e — > is characterized by 
that of the Wigner transform We{L, x, k, k). 
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3 The main results and assumptions on the random medium 



3.1 The main results 

This section presents our two main results. The first one describes the self-averaging of the back- 
propagated signal ijj^ in (|2^). Its proof is based on the second theorem, of independent interest, 
which shows that the Wigner transform converges in probability to a unique deterministic limit, 
solution of a transport equation. 

Theorem 3.1 Let the array function x(y) o.'^d the radially symmetric filter /(y) be in n L'^ {M.'^) , 
while Tpo G L^(]R'^) for a given k E M. Assume also that the refraction index n(z,x) satisfies 
assumptions in section \3.$ below. Then for each ^ G R'^ the back-propagated signal ip^ k) 



given by (2^) converges in probability and weakly in L^^(M'*) as e ^ to the deterministic function 



V^«(^,^;xo)= / e^'^<^-y^WiL,^o,^,K)My,^)^- (23) 



The function W satisfies the transport equation 

dW 1 

— + -k • VxW = kCW, (24) 

OZ K 

with initial data Wq{x,]s.) = /(k)|x(x)p and where the operator L is defined by 

f 7^^(^^^^,P-k)(A(p)-A(k)). (25) 



Here R{u!,p) is the Fourier transform of the correlation function ofV, defined by (Si.) below. 
The convergence of the Wigner transform is described by the following theorem. 



Theorem 3.2 Under the assumptions of Theorem 3.1 the Wigner distribution converges in 

probability and weakly in L^(M^°') to the solution W of the transport equation More precisely, 

for any test function A G L^(]R^'^) the process {We{z),X) converges to {W{z),X) in probability as 
e — > 0, uniformly on finite intervals < z < L. 



Here, (•,•) is the usual scalar product in L'^(M^'^). Theorem 3.1 is then a corollary of Theorem 3.2 



and Lemma |4.2| below. 

Theorem |3.1| provides a mathematical explanation for the important properties of the time- 
reversal experiment. First, it ensures that the back-propagated signal is essentially independent of 
the realization of the random medium in the high-frequency limit since k; xq) is deterministic. 

Second, it provides a quantitative description of the back-propagated field, which may be written as 

'i/'^(-,K;xo) = F(-,k;xo)*V'o(-,«^)- 

The kernel F{^, k; xq) is the Fourier transform in k ^ ^ of the solution of (^). We deduce from this 
formula that the back-propagated signal has much better refocusing properties in random media than 
in homogeneous media. As is explained in Q , the reason is that the limiting Wigner measure W is 
more regular in its k variable when the right-hand side in ( p4[ ) is present. This leads to better spatial 
decay of F and hence ipB localized tighter in a random medium than in a homogeneous medium. 
The improvement of the refocusing in random media has been carefully analyzed theoretically and 
numerically in |9|, in the Fokker-Plank approximation of (|2^). 
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3.2 Outline of the proof 



We summarize here the main steps of the derivation of the results stated in Theorems |3.1| and | 
We first obtain in Lemma 4.1 uniform bounds in for the Wigner transform We independently of 
the realization of the random medium. Lemma 4^ then shows that Theorem is a consequence 
of Theorem |3.2| . 

The main mathematical restriction is to assume that the random potential V{z, x) is Markovian 
in its first variable. We refer to section 3.3 for the details of its construction. Following our previous 
work in we show the convergence of the expectation E {(M^e, A)} — > {W , A) for every test function 
A. Its proof is simpler (and less general) than in because of the available a priori Lp' bounds on 
the Wigner transform, which satisfies the following Cauchy problem 

^ + k • V^W, = CeWe (26) 

M^,(0,x,k) =T^O(x,k), 



with 



{2-kY 



ty,(x,k-|)-VF,(x,k + |) 



Since k is a fixed parameter, we set k = 1 without loss of generality. The solution of (^6|) with initial 
data in is understood in the sense that for every smooth test function A(z,x, k), we have 

{We{z), \{z)) - {Wl A(0)) = f\We{s), + k • Vx + c)j \{s))ds. 

Here, we have used that is a self-adjoint operator for (•,•). Therefore for a smooth function 
Ao(x, k) we obtain (We(2:), Aq) = A£(0)), where Ae(s) is the solution of the backward problem 

d\ 

— ^ + k- VxA, + £eA,(s) = 0, 0<s<z 

OS 

with the terminal condition A£(2;,x, k) = Ao(x, k). This defines the process We{z) in L^(M^'^) 
and generates a corresponding measure on the space C([0, L]; L^(]R^'^)) of continuous func- 
tions in time with values in L?'. The measure is actually supported on paths inside a ball 
X = \W ^ LP' : \\W\\i^2 < with the constant C as in Lemma 4.L The set X is the state space 
for the random process We{z). We have proved in Q and will use in the sequel that the family P^ 
is tight: 

Lemma 3.3 The family of measures Pg is weakly compact. 

The proof of convergence of Ws to its deterministic limit is obtained in two steps. Let us fix a 
deterministic test function X{z, x, k). We use the Markovian property of the random field V{z, x) in 
z to construct a first functional Gx: C([0,L];X)— > C[0,L] by 

/•z ax 

Gx[W]{z) = {W,X){z)- / (VF +k.VxA + £A)(C)dC (27) 

Jo oz 

and show that it is an approximate P^-martingale. More precisely, we show that 

|E^^ {GxW\{^)\^s\ - Gx\W\{s)\ < Cx,lV^ (28) 
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uniformly for all W G C{[0, L]; X) and < s < z < L. Lemma |3.3| implies that there exists a 
subsequence £j — > so that P^. converges weakly to a measure P supported on C{[0, L]; X). Weak 
convergence of and the strong convergence (28) together imply that Ga[VF](^) is a P-martingale 



so that 

E''{Gx[W]iz)\J^s}-Gx[W]{s) = 0. (29) 



Taking s = above we obtain as in |Q] the transport equation (|2^ ) for W = {M^(-2)} in its weak 
formulation. 

The second step is to show that for every test function A(z, x, k) the new functional 

G2,x[W]iz) = {W,Xfiz)-2 / {W,X)iO{W, — + k-V^X + C\)iC)dC 

Jo 

is also an approximate Pg-martingale. We then obtain that {{W, A)^} — > {W, A)^, which implies 
convergence in probability. It follows that the limit measure P is unique and deterministic, and that 
the whole sequence Pe converges. 

That G2,a[W^](-2^) is an approximate P^-martingale uses very explicitly the uniform a priori 
bound on the Wigner distribution Wg. When the a priori bound is available only in a much larger 
space (the space A' of Q), we are not able to prove that the functional G2,a is an approximate 
Pe-martingale, and actually suspect that the result is not true. We expect to observe that some 
long-range correlations survive for sufficiently singular initial data, which would imply that the 
limiting measure P is no longer deterministic. 



3.3 The random refraction index 

We describe here the construction of the random potential V(z,x). The refraction index n(z,x) in 
non-dimensional variables is a random function of the form 

n2(z,x) = 1 - 2o-y(z,x), 

where the non-dimensional parameter a measures the strength of fluctuations. 

The mathematical analysis of the Wigner distribution in random media can be obtained by 
the method of diagrammatic expansions of the solution of the Schrodinger equation as in p^ . 
However, here we assume that the random field V{z, x) is a Markov process in z, as in Q, and we are 
able to analyze the Wigner transform directly. The Markovian hypothesis is crucial to simplify 
the mathematical analysis because it allows us to treat the process z i— > {V{z/e,x/e),Ws{z,x,'k)) 
as jointly Markov. 

In addition, V{z,x.) is assumed to be stationary in x and z, mean zero, and is constructed in 
the Fourier space as follows. Let V be the set of measures of bounded total variation with support 
inside a ball Bl = {|p| < L} 

V={v: [ \dV\ < C, suppVc Bl, Vip) = V*{-p)] (30) 

and let V{z) be a mean-zero Markov process on V with generator Q. The random potential V{z, x) 
is given by 

V(-.x)=/ 



9 



and is real and uniformly bounded: 

|y(z,x)| <c. 

We assume that the process V{z, x) is stationary in z and x with correlation function R{z, x) 

E {Vis, y)y{z + s, X + y)} = R{z, x) for all x, y G M*^, and z,s£R. 

In terms of the process V{z,p) this means that given any two bounded continuous functions (j){p) 
and '0(p) have 

E{{V{s)J){V{z + s),i;)} = {27ry [ dpR{z,p)(^{p)i^{-p). (31) 

Here (•, •) is the usual duality product on x M'^, and the power spectrum R is the Fourier transform 
of R{z, x) in x: 

R{z,p) = [ dxe-'P''R{z,x). 

We assume that R{z, p) G 5(M x M'') for simplicity and define R{uj, p) as 

R{uj,p)= [ dze~''^'R{z,p), (32) 
Jr 

which is the space-time Fourier transform of R. 

We assume that the generator Q is a bounded operator on L°°{V) with a unique invariant measure 
7r{V) 

Q*TT = 0. 

and that there exists a > such that if {g, vr) = then 

V'^gh^ < C|b||L-e-"^ (33) 

The simplest example of a generator with gap in the spectrum and invariant measure tt is a jump 
process on V where 

QgiV) = [ giVi)d7T{Vi) - g{V), [ d7T{V) = 1. 
Jv JV 

Given (^), the Fredholm alternative holds for the Poisson equation 

Qf = g, 

provided that g satisfies (vr, g) = 0. It has a unique solution / with (vr, /) = and ||/||l^ ^ . 
The solution / is given explicitly by 

f{V) = - / dre'^igiV), 
Jo 

and the integral converges absolutely because of (|3^). 

The particular Markovian model adopted in this paper is somewhat restrictive. However, the 
only information about the process that enters into the main result, the transport equation for the 
Wigner distribution in the limit e — > 0, is the two-point correlation function of V. Therefore, we 
expect that self-averaging of the Wigner transform and convergence to the solution of the transport 
equation hold for much more general classes of random perturbations of the refraction index. 
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4 Convergence of the Wigner transform 

4.1 A priori bounds on the Wigner transform 

The Wigner transform satisfies the following uniform bound. 

Lemma 4.1 There exists a constant C that is independent of e so that 



\We{z)\\L2(^^2i) < C||/||^2(Kd)||x 



(34) 



for all 2: > 0. 



Proof. It is shown in |21] that equation (^) preserves the L^-norm of W^: ||VFe(z)||i2 = ||We(2; 
0)11^2. The initial conditions for We are given by 



so that 



dydq 

(27r)'^' 



,ik-y+iq-y 



/(q)x(x- — )x(x + — 



dxdk\We{0,x,k) 



dydyidxdke 



ik-yi-ik-y 



/(y)/(yiMx-f )x(x + ^) 



xx(x-^Mx + ^) = (2^)'^ 



ey^ 



dxdy|/(y)|^ X(x-^ 



X(x+ yj 



and (34) follows. 

The above calculation also shows that the initial Wigner distribution W^(x, k) converges to 

l^0(x,k)=/(k)|x(x)|2. 

Notice however that one cannot expect strong convergence in L2 of Weiz) to W{z), the solution of 
the transport equation (24), because the latter has an L^-norm that decreases as z increases while 
the L^-norm of We is preserved. 

The next Lemma shows that one may drop the term e(y + ^)/2 in the argument of We in 
expression (^) for the back-propagated signal. 

Lemma 4.2 Let 0(xo) G L'^{R'^), then 



dxoV'f (^;xo)(/)*(xo) 



dxodke^^-^VF,(L,xo,k)0*(xo)^o(k) 



as e ^ 



(35) 



for all ^. 
Proof. We have 



V^f (^;xo)0*(xo)c?xo 



e*'^-^VF,(L,xo,k)0(xo)V^o(k)dxodk 



/ r (q)VF.(L, q,y - O^o(y) (e-'i (y+^)/2 - 1 



dyd<i 

{2'kY 



while 



/ 



l'^(q)nV'o(y)P 



1 - e 



dydq 



by the Lebesgue dominated convergence theorem. Therefore Lemma L2 follows from Lemma 4.1 . 

Theorem 3.1 then follows from Lemma 4.2 and Theorem 3.2, which is proved in the following 
sections. 
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4.2 Convergence of the expectation 



To obtain the approximate martingale property (28), one has to consider the conditional expectation 
of functionals F{W,V) with respect to the probability measure Pe on the space C([0, L];V x X) 
generated by V{z/e) and the Cauchy problem (^). The only functions we need to consider are 
actually of the form F{W,V) = {W, X{V)) with A G L°^{V;C^{[0, L];S{R'^'^))). Given a function 
F{W, V) let us define the conditional expectation 



T > Z. 



The weak form of the infinitesimal generator of the Markov process generated by Ps is given by 



dh w,v,z 



{W,X{V))]{z + h) 



1 



h=0 



d 



1 



X 



-{W,QX) + (W,(— + k-V^ + ^IC[V,-]]\ 



dz 



hence 



Gl = {W,\{V)){z) 



VF, ( + ^ + k • Vx + j] ) A ) {s)ds 



is a Pg-martingale. The operator K, is defined by 



/C[y,r,]V(x,77,k,y) = - 



i V (2vr)'^ 

The generator ( |36|) comes from equation (^) written in the form 

dW^ 1 

^ + k • VxVF, = ^/C[y(z/e), x/e]t^,. 

Given a test function A(z,x, k) G C^([0, we construct a function 

\e{z, X, k, y) = \{z, x, k) + VeAf (z, x, k, V) + eA|(z, x, k, V) 



(36) 



(37) 



(38) 



(39) 



(40) 



with Af 2(0 bounded in L°°(V; L2(]R2d)) uniformly in z G [0,L]. The functions Af 2 will be chosen 
so that 



\GlSz)-Gx{z)\\L^^v)<CxV~e 



(41) 



for all z G [0, L\. Here is defined by (^) with A replaced by A^, and G\ is defined by (|2^). The 
approximate martingale property (|2^) follows from this. 

The functions Af and Al are as follows. Let Ai(2;, x, r/, k, F) be the mean-zero solution of the 
Poisson equation 



k • V^Ai + QAi = -/CA. 



It is given explicitly by 
Ai(2;,x,r/,k,y) = - 



dre^ 



dV{v) 



3ir(k-p)+i{Tj-p) 



^(2)X,k- |) - A(2;,x,k + |; 



(42) 



(43) 



Then we let Af (z,x,k,l/) = Ai(z, x, x/e, k, y). Furthermore, the second order corrector is given by 
A2(z, X, k, V) = A2(z, X, x/e, k, y) where A2(2, x, 77, k, y) is the mean-zero solution of 



k-V^A2 + QA2 = /:A-Oi, 



(44) 
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which exists because E{/CAi} = £A, and is given by 



dre^ 



CX{z, X, k) - [/CAi] {z, X, 77 + rk, k, V) 



Using (H) and we have 
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h=0 



w^, ( ^ + k • Vx + 7] + 7^ ) + + ^^2) 



dz 



d 



WA — + k-V^jX + CXj + {W,i^^ + -k- Vx j (V^Af + eA|) + V^JC[V, ^]A| 
W^, ( ^ + k • Vx] A + £a\ + C^) 



with 



^ + k • Vx") Af + ("^ + k • Vx") A| + IC[V, ^]A|. 



The terms k- VxAf 2 above are understood as differentiation with respect to the slow variable x only, 
and not with respect to r/ = x/e. It follows that G\ is given by 



GUz) = {Wiz),X,) 



[ + ■ + ^) ^) " "^'^^^ ^^^^""^ ^^^^ 



and is a martingale with respect to the measure Pe defined on Z)([0, x V), the space of right- 

continuous paths with left-side limits [H]. The estimate (]28|) follows from the following two lemmas. 



Lemma 4.3 Let X £ C^{[0, L];S(M?'^)). Then th ere exists a constant Cx > independent of z € 
[0,L] so that the correctors Af(z) and X2{z) satisfy the uniform bounds 



and 



dXliz) 



dz 



+ k- VxAf(z 



^i(^)IIl°°(V;L2) + II-^2(^)IIl°°(V;L2) < Cx 

dXliz) 



+ 



dz 



+ k- VxA|(z) 



L°°(V;L2) 



<Cx. 



(46) 



(47) 



Lemma 4.4 There exists a constant Cx such that 

for any V £V and all e G (0, 1]. 

Indeed, (||) implies that \{W,X) - {W,Xe)\ < for all W e X and V e V, while 



Lemma 4.4 imply that for all z G [0, L] 



\\CHz)h.<C 



and 



(48) 



for allV eV so that (p|) follows 
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Proof of Lemma 4A. Lemma L4 follows immediately from the definition of /C, the bound ( |30D 
and the Cauchy-Schwartz inequality. 

We now prove Lemma f4.3| . We will omit the z-dependence of the test function A to simplify the 
notation. 

Proof of Lemma 4.3. We only prove (^). Since A G there exists a constant Cx so 

that 



|A(x,k)| < 



(1 + |x|5'i)(l + |k| 



5d-] 



Then we obtain using (30) and (|3 



\XUz,^,k,V)\=C 



\{z, X, k - |) - \{z, X, k + |; 







V 



<C I (ire-°^- sup/ |dl/(p)| [|A(z,x,k-|)| + |A(z,x,k+|; 



< 



C 



(l + |x|5'^)(l + (|k|-L)5'^X|k|>5L(k)) 



and the L^-bound on Ai follows. 

We show next that A^ is uniformly bounded. We have 



A|(x,k,y) 







Ai(x, ^ + rk, k - |, y) - Ai(x, ^ + rk, k + |, V) 



The second term above may be written as 



1 



,ip-{x/£+rk) 



Ai(x, ^ + rk, k - |, y) - Ai(x, ^ + rk, k + |, V) 



^^(P) Jp.(x/£+rk) 

(27r)^ 



^js{k-p/2)-q+j(x/£+rk)-q 



A(x,k----)-A(x,k-- + -, 



+ 



,ip-(x/e+r-k) 



dse" 



{2'kY 



,is(k+p/2)-q+i(x/£+rk)-q 



A(x, k + 2 - 2 ) ~ 22' 



Therefore we obtain 



|A|(x,k,y)| < C / dre 



V 



£A(x,k)| +sup / \dV{Y>)\ dse-'^'sup |(iVi(q 



Vi 



x(|A(x,k-|-|)| + |A(x,k-| + |)| + |A(x,k + |-|)| + A(x,k + | + | 



< C 



|£A(x,k)| + 



1 



(l + |x|5'^)(l + (|k|-L)5'^X|k|>5L(k)) 



and the L^-bound on Af in ( |4^ ) follows because the operator C : ^ is bounded. The proof of 
is very similar and is omitted. 
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Lemma 4^ and Lemma L4 together with (^) imply the bound ( p|) . The tightness of measures 
Pe given by Lemma |3l3| imphes then that the expectation E {14^£(2;, x, k)} converges weakly in L'^' 
to the solution W{z, x, k) of the transport equation for each z £ [0, L]. 



4.3 Convergence in probability 

We now prove that for any test function A the second moment E|(VFe,A)^} converges to (Vl^, A)^. 
This will imply the convergence in probability claimed in Theorem |3.2| . The proof is similar to 
that for E{(We,A)} and is based on constructing an appropriate approximate martingale for the 
functional {W iSi W, fj.) , where /^(z, xi, ki, X2, k2) is a test function, and W (Si W{z,xi,'ki,X2,'k2) = 
W{z,xi,\ii)W{z,X2,]<.2)- We need to consider the action of the infinitesimal generator on functions 
of W and V of the form 



F{W,V) = (M^(xi,ki)VF(x2,k2),Ai(^,xi,ki,X2,k2,T^)) = {W0W,fiiV)) 
where is a given function. The infinitesimal generator acts on such functions as 



h=0 



where 



2 1 

.7=1 ^ 



with 



and 



I 



dy(p)e'(p-^i) 



^{w^w,QX) + {w^w,n^2f^), 



/^(ki - |,k2) -/x(ki + |,k2) 



lC2[V,'n2]ii = - I dnp)e^(P-^2) ^(ki,k2-^)-Mki,k2 + ^) 



Therefore the functional 



G 



{w ®w,^l{v)){z) 



1 



d 



X2i 



(49) 



(50) 



(51) 



t^, ( -Q + ^ + ki • + k2 • + ^(/Ci[V^, — ] + K,2[V, {s)ds 

E OZ \l E E E 



is a martingale. We let /^(z, X, K) G ^(M^^ x M^f^^ ]-,g ^gg^ function independent of Y ^ where 
X = (xi,X2), and K = (ki,k2). We define an approximation 



iXe{z, X, K) = ix{z, X, K) + X, X/e, K) + e//2(X, X/e, K). 
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We will use the notation //f(2;,X,K) = X, X/e, K) and fj.2{z,'K,'K) = fi2{z,X.,'K/e,'K). The 
functions /ii and //2 are to be determined. We now use (p9|) to get 



' ■ dh 
1 



^^E^^y/iW ® M^, M.(V^))(^ + h) = -^(w®W,(^Q + J2y ■ V^. j /x^ (52) 



The above expression is evaluated at rjj = xj/e. The term of order e ^ in vanishes since n is 
independent of V and the fast variable ry. We cancel the term of order the same way as 

before by defining /ii as the unique mean-zero (in the variables V and rj = ir]i,r]2)) solution of 



It is given explicitly by 



^i(X,r7,K,T/) = - / dre^-« / dy(p)e"'(''i P)+*(^i-P) 



^ = 0. 



^(ki - |, ka) - ^(ki + |, k2) 



(53) 



+ - / dre'-'^ / dl/(p)e^^('^2-p)+i(r,2 P) L(ki, ks - ^) - /x(ki, k2 + ^) 
I . n .had L 2 2 



When /i has the form /i = A (8) A, then ^ui has the form ^ui = Ai (8> A + A (8> Ai with the corrector Ai 
given by (^). Let us also define ^2 as the mean zero with respect to vry solution of 



(54) 



(55) 



{Q + Y. k^' • v^.>2 + ^1 = E 1'^' 

i=i j=i i=i 

where / = / diryf. The function /_f2 is given by 

/>oo 

^2(X,r7,K,y) = - / dre^«[/Ci[F,r7i + rki]/ii(X,r? + rK,K) 

JO 

- [/Ci [y, + rki]^n] (X, 77 + rK, K, V)] 

- / dre^«[/C2[y,k2 + rT72]w(X,r/ + rK,K) 

JO 

- [/C2 V2 + ^k2]/ii] (X, T7 + rK, K, V)] . 

Unlike the first corrector ^1, the second corrector 112 may not be written as an explicit sum of tensor 
products even if ^ has the form ^ = \ ® \ because //i depends on V . 
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The P^-martingale G^f is given by 



G 



2,e 



{w (^w,^l{v)){z) 



Jo 



{W ®W,Ce)is)ds, 



where Ce is given by 



VF, ( ^ + ki • + k2 • + CD^i ) {s)ds (56) 



and the operator £1 is defined by 
1 



/ dpR{r,p) e'''(''i+ 2 ) P(/i(z, xi, ki, xa, ks) - /^(z, xi, ki + p, X2, kg)) 



2^" (27r)'^ 

e"-(ki-f)-P(^(^^xi,ki - p,X2,k2) - ^(z,xi,ki,X2,k2)) 



+ 



^^^l'2-P(M^,Xi,ki + E X2,k2 



/u(2;,xi,ki + f,X2,k2 + f )) 



+ 



^P-'^+'"''^-P(M2,Xi,ki - f,X2,k2 - f) - M^,Xi,ki - E,X2,k2 + E)) 
2*"^^-P+*'^-P(;U(^,Xi,ki - E X2,k2 + E) - ^(z,Xi,ki - E X2,k2 - f)) 

-e^''''^-P+^'^-P(M^,xi,ki + f,X2,k2 + E) - M^,xi,ki + E,x2,k2 - E)) 

e*''^''''+2)-P(^(z,xi,ki,X2,k2) - ^f(z,xi,ki,X2,k2 + p)) 
-e"'('^2"2)'P(/u(z,xi,ki,X2,k2 - p) - ;u(2;,xi,ki,X2,k2)) 



(57) 



We have used in the calculation of £1 that for a sufficiently regular function /, we have 



E 



dr e 



rQ 



dr 



R{r,p)f{r,p,-p)dp. 



The bound on (^s is similar to that on obtained previously as the correctors fij satisfy the same 
kind of estimates as the correctors A^: 

Lemma 4.5 There exists a constant > so that the functions fif 2 ofeey the uniform bounds 

II/^i(^)IIl2(R2'«) + llM2llL2(R2d) < Cf, (58) 



and 



dn\{z) 



dz 



+ ^k,-Vx,/iK^ 



L2 



5/i|(z) 



dz 



+ 5^k,--Vx,^|(z) 



L2 



<C„ 



(59) 



for all z e [0, L] and V GV. 



The proof of this lemma is very similar to that of Lemma L3 and is therefore omitted. 

Unlike the first moment case, the averaged operator £1 still depends on e. We therefore do not 



have strong convergence of the P^-martingale G^f to its limit yet. However, the a priori bound on 
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We in allows us to characterize the limit of G^f and show strong convergence. This is shown as 



follows. The first and last terms in (|57| ) that are independent of e give the contribution: 

dp 



poo p 

/ dr 

Jo JW' 



{2itY 



i2(r,p - ki)e*'' 2 (;u(z, xi, p, xg, kg) - ^(z, xi, ki, xg, kg)) 



, 2 2 



+i?(r, ki - p)e*'' 2 (^(z,xi,pi,X2,k2) - /i(z,xi,ki,X2,k2)) 

2 , 2 

+i?(z,p - k2)e*'' 2 (^(z,xi,ki,X2,p) - ^(z,xi,ki,X2,k2)) 
1 2 2 1 
k2 - p)e*'' 2 (/i(z,xi,ki,X2,p) - /i(z,xi,ki,X2,k2)) 



— ^-^,P - ki)(/u(2;,xi,p,X2,k2) - /i(z,xi,ki,X2,k2)) 



(27r 



+R{ ^-^'P - k2)(^(z,xi,ki,X2,p) - ;u(z,xi,ki,X2,k2)). 

The two remaining terms give a contribution that tends to as e — > for sufficiently smooth test 
functions. They are given by 



{27rY 



dr / dpR{r, p) x 



^,p.^2_2n+,,k2.p + girki.p+^^.p\ (^(^,xi,ki + f ,X2,k2 + f ) - /i(z,Xi,ki + f ,X2,k2 - f) 



We have 



i?(z,p) = i?(-z,-p)>0 



by Bochner's theorem. Since (£1 ~ ^2) and A are real quantities, we can take the real part of the 
above term and, after the change of variables r — > — r and p — > — p, obtain 



m - ^2)/u 



1 



(2vr) 



dr / dp ii(r, p) cos (p 



X2 - Xi 



-j|^gjrk2-p _|_ girki p-j 



(27r)'^ 



X (^//(z,xi,ki + E,X2,k2 + f) +/i(z,xi,ki - E,X2,k2 - f 
- ^i(z,xi,ki + f ,X2,k2 - f ) - ^(2;,xi,ki - f ,X2,k2 + f ) 

(ip(ii(-ki • p, p) + i?(-k2 • p, p)) COs(p • ^) 

X (^^(z,Xi,ki + f ,X2,k2 + f ) - /i(z,Xi,ki - f ,X2,k2 + f ) 
= 51 + 52 + 53 + 54 + c.c. 

We have (since /i is real-valued) 

I = (ixi(iki(ix2(ik2|(7i(z,xi,ki,X2,k2)|^ = C' / dxi(iki(ix2(ik2dp(iq-R(— ki • p, p)^(— ki • q, q 



xe 



*(P 'l)- ^a: ^//(z,Xi,ki - ^,X2,k2 + ^)/i(2,Xi,ki - ■^,X2,k2 + ■^). 



p 

2- 



q 

2' 



2' 



Using density arguments we may assume that fi has the form 

/i(xi,ki,X2,k2) = ^l(xi - X2)^2(xi + X2)/i3(kl)/i4(k2). 
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Then we have 



I = C (ixidki(ix2(ik2(ip(iqi?(— ki • p, p)/?(— ki • q, q) 

Xe-(P-'l)-^^^2(^^)^2(^^)^^(ki - |)^^4(k2 + |)/i3(kl - |);U4(k2 + |) 



e 



1/2 



= C||/i2|li2 / (iki(ik2(ip(iq^(-ki • p, p)^(-ki • q, q) 

X/i3(ki - -)^4(k2 + -)/^3(ki - -)/i4(k2 + -) 

where i^(x) = /ix(x). We introduce G(p) = sup;^ iZ(L(j, p) and use the Cauchy- Schwartz inequahty in 
ki and k2: 

|/| <C||/i2||i2||^*3|li2||/X4||i2 / dp(iqG(p)G(q) H^^) 
We use again the Cauchy-Schwartz inequahty, now in p, to get 

|/| <C||/i2||i2||/i3lli2||/i4|li2||G||i2 / dqG(q)f/ dpui?-)^) 

<Ce'^/2||/i2||i2||/i3||i2||^4|li2||G||,.2||G|Ui||z.||i2. 
This proves that IK^I ~ ^2)/^||l2 — > as e — > 0. Notice that osciUatory integrals of the form 

/ e*^^(p)c?p 

are not smah in the bigger space A', which is natural in the context of Wigner transforms and was 
used in [Q]. In this bigger space, we cannot control (£1 — ^2)fJ- and actually suspect that the limit 
measure P may no longer be deterministic. 
We therefore deduce that 

Gl = {W^ W, fiiV)){z) - j^' (w(g) W, + ki • V,, + k2 • + £2)^^ {s)ds 
is an approximate Ps martingale. The limit of the second moment 

W2{Z, XI, ki, X2, k2) = {W{Z, ^iM)W{z, X2, k2)} 

thus satisfies (weakly) the transport equation 

^ + (ki • Vx, + k2 • VxJW^2 = C2W2, 

with initial data VF2(0, xi, ki, X2, k2) = Wo(xi, ki)Wo(x2, k2). Moreover, the operator £2 acting on 
a tensor product A (8> A has the form 

C2[X(E> X] = CX(g) X + X(g) jC\. 

This implies that 

{W{Z, XI, ki)W^(z, X2, k2)} = E^ {W{Z, XI, ki)} E^ {W{z, X2, k2)} 

by uniqueness of the solution to the above transport equation with initial conditions given by 
Wq{xi, \<.i)Wq{x2,\<.2) ■ This proves that the limiting measure P is deterministic and unique (because 
characterized by the transport equation) and that the sequence VF£(z,x, k) converges in probability 
to VF(z,x,k). 
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